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Abstract 

DNA stretching experiments are usually interpreted using the worm-like 
chain model; the persistence length A appearing in the model is then in- 
terpreted as the elastic stiffness of the double helix. In fact the persistence 
length obtained by this method is a combination of bend stiffness and intrinsic 
bend effects reflecting sequence information, just as at zero stretching force. 
This observation resolves the discrepancy between the value of A measured in 
these experiments and the larger "dynamic persistence length" measured by 
other means. On the other hand, the twist persistence length deduced from 
torsionally-constrained stretching experiments suffers no such correction. Our 
calculation is very simple and analytic; it applies to DNA and other polymers 
with weak intrinsic disorder. 
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Introduction and Summary: The DNA in living cells is often described as a passive database 
of pure information, the genome. In fact, however, the DNA molecule itself actively col- 
laborates in its own packaging, transcription, regulation, and repair Unraveling the 
underlying mechanisms of these crucial processes requires an understanding of the basic 
mechanical properties of the DNA duplex. For example, the fundamental unit of DNA 
packaging, the nucleosome, is delicately balanced between elastic stresses and bonding ener- 
gies 0; an accurate account of the former is clearly important for analyzing the stability of 
the whole complex. Since nucleosomal DNA is under torsional as well as bending stress |§ , 
an accurate model incorporating both twist and bend is needed. 

Recently a new class of experiments has permitted precise physical control over single 
molecules of DNA [||. For example, a single molecule of known contour length L can be 
subjected to known stretching force / at its ends and the resulting extension (end-to-end 
length) Z measured. Simple arguments from polymer physics then predict that Z < L since 
thermal fluctuations keep a flexible rod from being perfectly straight; Z approaches L at large 
/. Remarkably, Bustamante et al. found that a very simple model, the "worm-like chain" , 
fit the force-extension data over four orders of magnitude in / ||. The model attributes 
to DNA just one parameter, the bend persistence length A e g; subsequent experiments have 
refined its value to || A e s = 40 nm [0. In a refinement of the technique, Strick et al. 
devised a tor sionally- constrained stretching experiment [§|§; analyses of the corresponding 
directed walk problem led to values of the twist persistence length C e g between 75 and 
120 nm |]10|-[T3[]; in each case C e s/A eS was found to exceed unity. 



The purpose of this note is to show that the value A eS measured by stretching experiments 
does not directly reflect the bend stiffness of the DNA helix, but rather a certain combination 
of stiffness and disorder induced by the sequence of natural DNA. Since these two effects will 
enter in different combinations in other circumstances, for example the nucleosome binding 
energy, it is important to disentangle them. In fact A c g underestimates the true elastic 
stiffness A, while C e s accurately reflects the true C, as announced in Thus the large 



observed value of C e s/A e g is perhaps not as mysterious as it at first seems. 



Recently Bensimon et al. have independently studied these and other issues |14[ . Using 
a different model from ours, they found analytical formulae for low-force stretching and 
numerical results for all /, at both strong and weak disorder. Below we will restrict to the 
case of weak disorder, the case relevant for DNA. In this limit the calculation becomes very 
simple. The result obtained here for A e g differs from [|14|| , as described below. We will also 
retain the torsional degree of freedom needed to study the twist stiffness. 

The result of this note is perhaps not surprising in the light of extensive earlier work 
on DNA coils at zero applied tension. A uniform rigid stack of monomers must form some 
sort of helix, and in particular such a helix will have a straight axis in its undeformed state. 
DNA, however, is a stack of four different types of unit. The sequence of natural DNA has a 



small component with period equal to the helix repeat |[15|| , but mainly the sequence imparts 
random natural bends to the rod [|TIJ . Trifonov et al. noted that even in the absence of any 
thermal fluctuations a randomly-kinked rod would follow a random walk of some persistence 
length P, which they called the "static persistence length." They argued that the effective 
persistence length of such a coil at nonzero temperature would be |T7| 



A cS = A/(l + \), (1) 

where A ■ k^T is the true elastic stiffness of the rod and A = A/P, and they verified 
formula ([I]) with Monte Carlo simulations fli~8" |. Trifonov et al. computed the numerical 



value P = 216 nm and hence A = 0.3 starting from sequence information and estimates of 
the wedge angles. Later Bednar et al. measured A more directly by comparing random coils 
of natural DNA to synthetic constructs designed to be straight; they obtained A = 78 nm, 
A eS = 45 nm, and hence A = 0.4 [K$ZQ |2|. 



One might imagine that under extensional force the kinked rod would simply follow 
the usual worm-like chain result with A replaced by A e g from (|1|). Indeed this is correct 
for weak disorder (small A). In contrast, Bensimon et al. found that at weak disorder 
A eS = A{1 - \y/X) @, while Marko and Siggia argued that at high force disorder is 
immaterial: A e g = A ||22||. 



For a rod under torsional stress similarly define C c g by the torque r needed to change 
the linking number Lk to a value different from its relaxed value Lk : 

Lk - Lk 



r = C cS ■ k B T ■ u - 



Lk 



(2) 



Here ujq = 1.85/nm is the rotation per unit length of relaxed DNA. In contrast to ([!]) we 
have 



C, 



cff 



c + - 



(3) 



where the ellipsis denotes terms vanishing at high force or greater than first order in A. 
Calculation: We wish to evaluate the extension of a randomly-kinked, flexible rod under 
an imposed tension /, and later an applied torque as well. We seek the leading term in an 



expansion in weak disorder; the extension to higher orders is straightforward [23 ). 

To describe the rod conformations, let E a (s) be an orthonormal triad describing the 
orientation of the rod segment at arclength s from the end, with E% the tangent to the rod 
axis. The spatial components E ia of these three vectors thus form an orthonormal matrix 
E(s). Let Q = E _1 E = J2i ^«Tj, where the dot denotes d/ds. Tj are the three antisymmetric 
3x3 matrices generating rotations, e.g. [Ti] 2 3 = +1. The elastic energy of a conformation 
is then: 



^eias/ k B T = - ds 
2 Jo 



(4) 



To this energy we now add a term describing the work done by the external force, 

/ 



dsE. 
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(5) 



The functions Ci( s ) appearing in (^) specify the random kinks p4 |. We give them an 
isotropic, Gaussian distribution: 



= 0; = j5(s - s') 



(6) 
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Here the double brackets signify an average over an ensemble of many possible sequences 



25| . Considering the curve whose curvature is exactly Ct( s ) one can see th & t P = A/X is the 
structural persistence length mentioned above, by calculating [^ 3 (0)-£ l 3(s)] = 1 — -p + 0(s 2 ). 
The constant g in eqn. (|6|) will drop out of our answers. 

Thus even neglecting thermal undulations altogether, straightening the rod requires some 
extensional force to overcome the elastic energy (|4]). We must now introduce entropic effects 
as well, and compute the full extension 

Z/L = [(E 33 (0))j , (7) 

where the angle brackets are the usual thermal average. 

To carry out the calculation, begin with the Euler angle representation of a rotation 
matrix, defining three fields 9(s), and ip(s) by 

E = g-^e-^e"^ 3 . (8) 

To exploit the assumed isotropy of the rod and its disorder, define the complex variable 
W = (fii + in 2 )/V2 = e-'^i-W + 0sin0)/V2. Similarly let Z = (& + ^Vv^, which then 
obeys [Z(s)Z*(s')J = -jS(s — s') and [Z(s)Z(s')J = 0. The energy then becomes 

./ 



E e i as /k B T = J ds 



A{\W\ 2 - WZ* - W*Z) + \C{i) + 0cos# + Ca) 2 - r^costf 

KB J- 



■ (9) 



We have dropped the divergent constant / \ Z\ 2 from (||) because constants in the energy do 
not affect thermal averages. It is now clear that the disorder field Cs( s ) may be eliminated 
from the last term of @ by shifting the definition of ip. Since ip does not enter the first 
term, while the next two terms already contain the disorder field Z, this shift eliminates £3 
altogether to leading order in the strength A. The physical meaning of this shift is simple. 
Consider a straight, isotropic rod with a randomly- rotating reference stripe painted on its 
surface. Nothing changes if we pass to a different reference frame rotated at s by an angle 
J s ds' (3(3') relative to the old one. 
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What makes our problem interesting is that the disorder Z can not be so trivially elim- 
inated, due to a clash between the A terms and the / term. To leading nontrivial order in 
the disorder strength A the JJ-terms of (|9|) contribute 

1 + A 2 J dsds' W(s)Z*(s)W*(s')Z(s') , (10) 

to the Boltzmann weight e~ Eclas ^ kBT . Performing the average over the Z fields eliminates 
one of the integrations over s, so that the correction factor is the leading term of e AX I ds l w l 2 . 
Comparing to (fj|), we see that to 0(A) the effect of disorder is simply to replace A by 
A e ff = A(l — A), leaving C unchanged. This proves (0,^) since we are working to first order 
in A. To go beyond this order we must be careful to treat the disorder as quenched, for 
example via the replica trick . 

We can easily incorporate an external torque r applied at the ends of the rod: r couples 
to the change in Link density, which in our variables is simply ip + <fi + We added (3 to 
the formula of |27| in order to measure the change in Link from the unstressed value; the 
same shift in the definition of ip used earlier thus eliminates £3 here as well. 

Thus within our approximations the only effect of sequence on entropic elasticity is to 
reduce the effective bend persistence length, as claimed in eqns. (0,0)- The first correction to 
eqn. (H) in powers of 1/ v7 is & lso simple to obtain by substituting eqn. ([!]) into the formula 
for the effective stiffness given in [p~2| , [3~3j| ||28|| . 

Discussion: The model investigated above may seem highly reductionist, neglecting as 
it does all the specific properties of DNA, e.g. the specific bends at particular base-pair 
junctions. Indeed we have used a continuum model, where there are no base-pairs at all. 
But it is precisely the existence of a good continuum limit, despite the very singular form 
of the assumed disorder (^), which gives the result universality. Like the phenomenon of 
entropic elasticity itself, random kinks affect the force-extension curve via fluctuations over 
length scales much longer than a base-pair. 

The analysis given here explains the qualitative success of models without disorder in 
fitting DNA stretching experiments. It also predicts that single-molecule stretching ex- 



periments on long, intrinsically-straight DNA would show the same increase in effective 



persistence length seen at zero force, for example in [|H]|2(J. More importantly, it implies 
that the elastic stiffness relevant for deformation of a given segment of DNA on scales shorter 
than a micron is considerably greater than the value obtained by fitting the worm-like chain 
model to stretching experiments. 
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